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Abstract 

The contributions to the coefficient functions of the quark and the mixed quark- 
gluon condensate to mesonic correlators are calculated for the ffist time to all 
orders in the quark masses, and to lowest order in the strong coupling constant. 
Existing results on the coefficient functions of the unit operator and the gluon 
condensate are reviewed. The proper factorization of short- and long-distance 
contributions in the operator product expansion is discussed in detail. It is found 
that to accomplish this task rigorously the operator product expansion has to 
be performed in terms of non-normal-ordered condensates. The resulting coef- 
ficient functions are improved with the help of the renormalization group. The 
scale invariant combination of dimension 5 operators, including mixing with the 
mass operator, which is needed for the renormalization group improvement, is 
calculated in the leading order. 
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1 Introduction 



Since the pioneering papers by Shifman, Vainshtein, and Zakharov (SVZ) |T[], 
QCD sum rules have played a major role in extracting parameters describing 
the structure of the QCD vacuum, the so called QCD condensates, as well as 
for calculating hadronic parameters, e.g. masses and decay constants of mesons 
and baryons from the fundamental QCD Lagrangian. This procedure assumes 
the "duality" of a hadronic versus a partonic description, the partons here being 
quarks and gluons. 

The calculation of the partonic part is performed in the framework of the 
operator product expansion (OPE) 0, most commonly including operators up 
to dimension 6, whose contributions signal the breakdown of perturbation theory 
for low momenta and parameterize non-perturbative effects. A large amount of 
results on the coefficient functions in the OPE exists in the literature, usually 
invoking some additional approximation, as for example small or large quark 
masses, or equal masses in the case of mesonic correlation functions. For further 
information the reader should consult refs. and 0-0. 

An important aspect which has to be mentioned in this context is the proper 
factorization of short- and long-distance contributions in the OPE. This means 
that the coefficient functions describing the short- distance part of a given corre- 
lation function should be free from dependences on the infrared structure of the 
theory, e.g. mass logarithms, and all such dependences should be absorbed into 
the corresponding matrix elements of operators, the condensates, which shall con- 
tain the long-distance contributions. This "cancellation of mass log's" has been 
already discussed in the literature [jlO|, 0, |12|. We shall deal further with this 
point below. 

In this work we present new results on the coefficient functions for the quark 
and the mixed quark-gluon condensate in the unequal mass case, to all orders in 
the quark masses, and to lowest order in perturbation theory, for scalar, pseu- 
doscalar, vector, and axialvector mesonic correlators. These results are obtained 
without any of the approximations mentioned above. In addition, for complete- 
ness we review the corresponding results for the unit operator and the gluon 
condensate. Extensive comparison of our results with results existing in the 
literature is made in various limiting cases. We explicitly show that all mass log- 
arithms, which for small masses correspond to long-distance contributions, can 
be absorbed into the QCD condensates only, if the OPE is performed in terms 
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of non-normal-ordered operators. This observation was already made in []T3|, ^^ . 



For an application see also ref. [T^. The resulting coefficient functions are im- 
proved with the help of the renormalization group. In this context, the scale 
invariant combination of dimension 5, including mixing with the mass operator, 
is calculated. 

An application of our results to an improved determination of the current 
strange quark mass, as well as a discussion of higher order as corrections, will be 
presented in a forthcoming publication [p!6| . 

In sect. 2, we summarize some known facts on the OPE. The new results 
for the coefficient functions of the quark and mixed condensate are presented in 
sects. 4 and 6, and the corresponding results for the unit operator and the gluon 
condensate are given in sects. 3 and 5. The proper factorization in the OPE and 
the cancellation of mass log's is discussed in sect. 7, and in sect. 8 we show how 
to improve the coefficient functions with the help of the renormalization group 
(RG). Finally, sect. 9 summarizes our results. 

2 The Coefficient Functions 

Let us first define the relevant two-point and coefficient functions which we are go- 
ing to examine, and then discuss some of their general properties in the framework 



of the OPE g, . The scalar, pseudoscalar, vector, and axialvector two-point 
functions are defined to be the matrix elements, between the physical vacuum 
of the time-ordered product of the corresponding currents, 

n^(g) = t{n\TN,{jriq)jlm\n), (2.1) 

where F stands for one of the Dirac-matrices F G {1,^75,7^,757^}, specifying 
the quantum numbers of the current (S, P, V, A respectively), the tilde sign 
will always denote quantities Fourier-transformed to momentum space, and Na 



is a suitable renormalization procedure for the operator product [T^, |T8|, |T^. 



Throughout this work we will use dimensional regularization in D = 4 — 2e 



dimensions [^, |2T| and the modified minimal subtraction (MS) scheme |52 
The current jr shall take the form 

Mx) =:Q{x)Tq{x):, (2.2) 

where Q and q are quarks of possibly different fiavour with masses M and m 
respectively. If M and m differ, we shall always assume M > m. To the or- 
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der we are working, the normal-ordering is sufficient to renormalize the current. 
However, in higher orders, additional subtractions are required. 

Using the Ward-identity relating the divergence of the vector (axialvector) 
two-point function to the scalar (pseudoscalar) two-point function |2^, we 
can express n^(g) in terms of 4 Lorentz-scalar two-point functions 11^ (g^), with 
/ G {S, P, V, A}. Then the vector (axialvector) two-point function can be written 
as 



(M =F m 
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+ 9,.^-^^[{QQrT{qqf^], (2.3) 

with Il^'^{q'^) being the scalar (pseudoscalar) two-point function. In the following, 
(...)^^-* will always denote matrix elements between the physical vacuum | f2) 
(condensates), the superscript indicating renormalization to the order considered 
(1-loop in this case). Please note that the quark condensates in eq. ( p.3|) are 
not normal-ordered. The reason for this will be thoroughly discussed below. On 
the other hand, because of this fact, the third term in eq. (|2.3|) depends on the 



renormalization scale — m{qq) is no longer renormalization group invariant |]T4 
— but this dependence is cancelled by the second term. 

In the framework of the OPE, T{^{q^) can be expanded in terms of local 
operators Oj times Wilson-coefficient functions Cj, Ci denoting the coefficient 
function of the unit operator, 

tf(g^) = C;(g^M,m,/.^a.) + E ^i^^^^^ (2-4) 

thereby separating the short- distance dynamics, described by the Wilson-coefficient 
functions, from the long-distance behaviour, included in the operators Oj. The 
bar indicates coefficient functions corresponding to tree-level matrix elements, ^ 



is a renormalization scale in the MS'-scheme, and the rii are the canonical di- 
mensions of the Oi. With [{rii — l)/2] denoting the integer part of {rii — l)/2, 
it follows that the c/ have dimension or 1, depending on the dimension of Oj 
being even or odd. 

Proper factorization of short- and long-distance contributions in the OPE re- 
quires the calculation of the matrix elements in eq. (|2.4|) at the same order as 
the coefficient functions. This procedure is well known from deep inelastic scat- 
tering [P^, El] and weak decays P31. In this way, taking into account mixing 
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of the operators under renormalization, all dependences on the infrared struc- 
ture of the theory are absorbed into the matrix elements. Calculating the C- 
straightforwardly using Wick's theorem, naturally leads to normal-ordered op- 
erators. However, as has been proven by p6| , p7| , |28[| , in minimal subtraction 
schemes the coefficient functions are analytic in the masses, if the OPE is per- 
formed in terms of non-normal-ordered operators. This feature is in general not 
present in non-minimal schemes. Since normal-ordering is a non-minimal scheme, 
using normal-ordered operators there appear terms of the form m** In^'m^///^ in 
the coefficient function Ci. These terms would make it impossible to calculate 
non-leading contributions, because, while summing up the usual perturbative log- 
arithms by setting //^ = — g^, we would acquire possibly large log's of the form 
In {w? / — q^). For small masses they are remnants of the long-distance structure 
of the theory and should be absorbed into the corresponding condensates. There- 
fore, to arrive at physically sensible coefficient functions, we have to express the 
tree-level matrix elements (: Oj :)^^^ in terms of the renormalized {Oi)'^^\ and, 
working with non-normal-ordered operators, also mixing with the unit operator 
under renormalization has to be taken into account. This will be done below. 
In our analysis we shall include operators up to dimension 5, namely 



{QQ). m. i-FF), (gsQaFQ), (gsqaFq) } , (2.5) 



related to the quark condensates, the gluon condensate and the mixed conden- 
sates having dimension 3, 4, and 5 respectively. Implicitly, summation over 
colour, spinor, and Lorentz indices is assumed, and the a-matrix is defined as 
<7fj.u = i/2 [jfj.y'yu]- Coefficient functions in various approximations up to dimen- 
sion 8 have been calculated, and can be found for example in refs. [Q, @-|^, 



and 12, 29, 30 



Calculating the diagrams of fig. 1, in the leading order one obtains the fol- 



lowing relations [jlO|, |T4| 



47^2 



M' 



2 



/i2 



1 



,(0) 



12M ' TT 



= {^:FF:f\ (2.6) 

vr TT 



{g.QaFQi^.)f^ = (:.g,QaFQ-f^ + ^ln^{^:FF-f^ 
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where is the number of colours. The corresponding relations for the g-quark 
condensates are obtained through the replacements Q ^ q and M ^ m. Since 
the gluon condensate is already 0{as), it does not get renormalized to the order 
we are working. 

Inserting these relations into eq. ( |2.4|) , we find the infrared stable coefficient 
functions C/, 

Ol - Cl.^{f\u.^-Acl,,.4\"^4-ACL], (2.7) 



q2 



^OQ - ^QQ ^Iq - *^<7a ' (2- 



qq ^ qq 



.2 



^FF - ^FF + + 12 m ^^"-^ 2g2 ^QFQ ^2 ^^Fq , 

(2.9) 

^QFQ = ^QFQ C|pg = Cg^g . (2.10) 

In terms of these coefficient functions, the OPE takes the form 

n^(g2) = Cl{q\ M, m, a,) + ^ ^^^f^;^:!!;)/,/ ■ (2.11) 

In the following sections we calculate the coefficient functions of eqs. (|2.7| )- 
( |2.10D for the unit operator and the quark, gluon, and mixed condensate to all 



orders in the quark masses, and show explicitly that up to operators of dimen- 
sion 7, they are analytic functions of the masses. Of course, as we shall also 
demonstrate, the coefficient functions may still contain non-analytic pieces of 
higher dimension, which are only cancelled through mixing with higher dimen- 
sional operators. In our example this is the case for Cf and Cpp. This point will 
be discussed further in section 7. 



3 The Perturbative Coefficient 

The perturbative coefficient function for two different quark masses has been 
already calculated in ref. [||, in the leading order, as well as to 0{as). For 
completeness, and for further reference, we give here our results for the coefficient 
functions c/, which are in agreement with however using a slightly more 
compact notation. 

Ci'^'ig') = -£-^[qlls{q',M,m,^i')-M'lM-m'l^Y (3.1) 
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N 
12^ 



2\2 



(M2 - m2) 



+ |. _ M2 - ^ . (3.2) 



Here Im = ln(M^//i^) — 1, and Isiq"^, M,m, is the scalar one-loop integral, 
Is{q\M,m,^?) = 



1 , , xM'^ + (1 — x)m'^ — x(l — x) q^ 
ax In 



, u+l (M'^-rn^), M , Mm 
In + ^ In h In — — 



with 



q% = - (M ± m) 



u-l 



q 



m 



and 



u = Jl - 



4Mm 



-2, (3.3) 



(3.4) 



To show explicitly the appearance of mass logarithms in the coefficient func- 
tions, we also present their expansions in the quark masses up to order m^. 



N 



2g|-(g|-M2-m2) In ■ 



-q 



- ± 2M^m + 2M2m2 ± 2Mm^ + -m^ 
2 2 



+ 



M m 
(M ± 2m) In — - + (m ± 2M) In 



2 ^ 1 



N 



127r2 I 3 



-9 



-r 



(3.5) 



r 1 1 A/f^ 

-M^-2MW + -m^ + M^ Ini^ + m^ In -3 3 ^. (3.6) 



2 
2 



1 



_g2 J q2 



As has been discussed in the previous section, these mass logarithms are cancelled 
after inclusion of the mixing with the quark condensate. This will be performed 
in section 7. Expressions for the c/, only expanding in m, are provided in 
appendix A. 



4 The Quark Condensate 

The calculation of the coefficient functions for the quark condensate follows 
closely the method used in ref. PT|. The contribution of the quark condensate 
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to the two-point function 11^ (q), shown in fig. 2, is given by 

K+QQi^) = - J 7^ (:?(0)r^Q(p - q)Tq{p) -f^ + {q ^ Q) . (4.1) 



Here Sq{p — q) = (^ — ^— M)~^ is the free quark propagator. 

A necessary ingredient for calculating the coefficient functions to all orders in 
the quark masses is a closed expression for the non-local quark condensate. In 
x-space this expression reads |3l|, ^ 

{■MO)q,ix):)f^ = (: qq :) T f (z ^ + m),. f tT^""^ ^f^^^,, , (4.2) 



4m''^^'' \2j''^ "'°^on!r(n + D/2) ' 

where the index qq denotes the projection onto the local quark condensate. One 
should remark that the non-local quark condensate has, of course, additional 
contributions from higher dimensional operators (see section 6). The sum in 
( [4.2|) can be expressed in terms of Bessel functions, 

~ (-mV/4)" / 2 \''^''\ ^ .^ 



but like in [^, we prefer to work explicitly with the expanded form. Since we 
shall derive similar relations for the mixed condensate, we skip the derivation of 



eq. ([4.2|), and refer the reader to appendix A of ref. p2 



The corresponding p-space expression for the non-local quark condensate is 
given by ||5T[ 

(4.4) 

It is easy to verify that the quark condensate ( [4.4|) satisfies a free equation of 
motion, 

(^-m)(:g(0)g»:)(y = 0. (4.5) 

This is very useful, as it allows to replace an arbitrary non-singular function 
/(j9^, j9, . . .) by /(m^,p, . . .) in integrals of the type 

'd''pfip',p,...){:qmip)--)f^, (4.6) 
thereby greatly simplifying the calculation. 
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After performing the momentum integration and using the relations 



m 



4 P 



{2n) \ {m?q 



(4.7) 



2Fi(l,l/2;D/2;4z) j- 



1 - VI -4z 



(4i 



where 2F1 is the Hypergeometric function |3^, , we arrive at the following final 
results for the two-point functions in Z) = 4 dimensions: 



2m 
3 m 



r 



A 

M2 + m2 



(4.9) 



1 + 2 



(M2 - m2) 



[ + + - 2 (M^ - m^)^ /g^ ] 
[g2 _M2 +m2] 



with 



2 2 
m q 



2 12 ■ 



(4.10) 



(4.11) 



The corresponding functions for the heavy quark Q, n^Q(g^), can be obtained 
through the replacements q ^ Q and m M. In the equal mass case, M = m, 
our results agree with the various results given by the authors of refs. P, ^ 
except for the axial current in 0. 

To examine the structure of the coefficient functions more explicitly, and to be 
able to compare with other previous results, we expand the two-point functions 
in the quark masses up to operators of dimension 6. For the coefficient functions 
of the quark condensate this leads to 



m 

^M--T — 
2 



(4.12) 



-.VA, 



m 



1 + 2 (m^ - 



3 'q 



(4.13) 



and related expressions for the Q-quark. Again, in appendix A we provide results 
for solely expanding in m up to order rn? . These results also agree with the results 
given by Generalis 0. 
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5 The Gluon Condensate 



The contribution of the gluon condensate to all orders in the quark masses has 
already been presented in refs. [|, 0]. We agree with these results and, for 
completeness, and further reference, cite here the corresponding expressions. 

^ (^■.FF-f\ (5.1) 



l2Mm ' TT 



48 ' TT ■ I 2n5 m - 1 

(5.2; 

For P and A an additional change M — >■ — M in m, being equivalent to u —>■ 1/u, 
has to be performed. In the equal mass case our results for Upp agree with those 
of Bagan et.al. 0. 

We would like to point out that the non-transverse part for the two-point 
function np]S|^^ of ref. |Q (eq. 3.53), and our scalar (pseudoscalar) two-point 
function Il^'j^, differ by the g-independent piece in eq. Writing eq. (|2.3|) in 

terms of tree-level condensates would cancel this additional term, but then the 
Ward-identity would no longer be valid. 

Upon expansion in the quark masses up to operators of dimension 6, we find 
for the coefficient functions Cpp 

-^SP 1 1 /M m\ 1 , r, r,, 

24 12 Vm M/ 12g2^ ^ 
1 /M^ Mmr , Mm 



3 + 21og^ , (5.3) 



12g2 Km MJ Aq 

Clf = + (5.4) 

Like in the case of the quark condensate, in appendix A we provide results for 
the Cpp, expanded only in m. 

A few remarks are in order here: as is evident, at this intermediate stage 
there appear 1/m as well as m^logm terms. Nevertheless, they are remnants 
of the long-distance structure of the vacuum condensates and will cancel in the 
final result for the coefficient functions Cpp, once the additional contributions 
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through mixing ( p.9| ) have been included. This cancellation of mass singularities 
has already been extensively discussed in |10, 11, 12]. We present the coefficients 
Cpp in section 7. 



6 The Mixed Condensate 

The calculation of the contribution to the mixed condensate to all orders in the 
masses is somewhat more complicated than in the case of the quark condensate. 
We shall therefore discuss its evaluation in slightly more detail. The two diagrams 
contributing to the mixed condensate are shown in fig. 3. 

Let us ffist consider diagram 3a. Working in the coordinate gauge, this graph 
also arises from eq. ( [4.1| ), since in the interacting case q{x) has to be expanded in 
terms of covariant derivatives (see for example p6|), and hence yields contribu- 
tions which involve gluon fields. The procedure for the calculation of the mixed 
condensate contribution is similar to the one presented in the appendices A and C 



of ref. [0. For the explicit steps, the reader is referred to this publication. Our 
result for the non-local quark condensate contribution to the mixed condensate 
is 



OO 

n=0 



(n-l)i^ + nm] \ ' { (6.1) 



Up to terms 0{m?), the order calculated in ref. our result agrees with 

Similar to the quark condensate, the mixed condensate contribution satisfies 
an equation of motion. In p-space, this equation of motion reads 



(p^ - m^) (: mm = - ^'-'S^'S ^"S ' (6-2) 



which also allows for simplification in momentum integrals. The replacement in 
momentum integrals which can be made by means of eq. (|6.2|) has the form 



rf^p/(p^p,...)(:g(0)g(p):)(°^ j d^p \ /(m^,^, . . .) (:g(0)g(p) :)(°^ 



■.gsqaFq-.r r5^,^(^2^^^___)l (:^-(o)g»:)(;) I. (6.3) 
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The contribution from diagram 3b can be evaluated by considering the inser- 
tion of one gluon field strength into the non-local quark condensate of eq. (^4.2| ) 
[ P^ . The corresponding expression for the non-local mixed condensate is found 
to be 

(:Ma(0)F,.(0)g,(x) = ^^j^ _ ^^^j^ _ ^2 i'- 9sqcrFq-f^ T 

ht^d, - 7,9^) + m 0-^, J (i ^ + m)^^ J2 n\T{n + D/2) " ^^'^"^ 

Here = t"-F'^j^, where are the generators of the colour group SU{N). The 
expansion of eq. ( |6^ ) has been calculated by the authors of ref. up to order 
m^, and we agree with their result, except for the term 0{m^). The derivation 
of eq. (|6.4|) is presented in appendix B. 

It is obvious from the structure of eq. ( |6.4| ), that the non-local mixed con- 
densate also satisfies the free equation of motion 

{^-m){:g,UO)F,Mqp{x)-)% = 0. (6.5) 

Using eqs. (|0| ) and (|6.4| ), as well as the corresponding equations of motion, 
we obtain the following final results for the contribution of the mixed condensate 
to the two-point function of eq. (|0|): 

{■.gsqaFq:)^^^ 



nri4(^^) = -%g^{[g^ + 2M(MTrn)](g^-M^)^-rnV(g^ + M^) 
o III q q+q^sz ^ 



f{z^) , (6.6) 



- 2Mm\M T m){2M- - m^) - } , (6.7) 

with 

P(g^ M, m) = + 2M(M t m)] (g^ _ - M^m'^q^{4:M^ t GMm + m^) 

- m^q\q'^ + M^) + 2Mm'^{M T m)(3M^ - 3M^m^ + m^) . (6.8) 

As in the case of the quark condensate, the corresponding functions for the heavy 
quark Q, UQp^^q'^), can be obtained through the replacements q ^ Q and m ^ 
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M. In the equal mass case, M = m, our results agree with the result given in ||^. 
The expansion in one light mass m up to operators of dimension 6 for the vector 
current agrees with the result obtained in [0]. 

Expanding the coefficient functions for the mixed condensate up to operators 
of dimension 6 yields the well known expressions 

C^fV) = ctfM") = ±f ' = C^^,%') = 0, (6.9) 

and related expressions for the Q-quark. Again, in appendix A, we give results 
for the expansion in one small mass m. 



7 The Cancellation of Mass Log's 



We are now in a position to calculate the infrared stable coefficient functions 
C[ and Cpp of eqs. (|2.7|) and ( ^l9|) . The result up to dimension 4 for the unit 
operator is found to be 



Cf'V) = #^(4g|-(M^ + 4MW + m^)l 



-2 



1 1 



qi- -m' + (M^ ± 2M^m ± 2Mm^ + m^) , 

g2 



In- 



(7.1) 



N 
72^ 



10 + 18 (M^ + m^) - 9 ( 3M^ - AM'm^ + 3m*) — 



-6 



g2-3(M* + m* 



In 



(7.2) 



and up to dimension 6 for the gluon condensate 



1 (M^ + m^) Mm 



1 

12 



12 g2 

(M^ + rri^) 
18g2 



± 



4g2 



3 - 21n 



(7.3) 



(7.4) 



It is evident that all mass logarithms have cancelled, and the resulting coeffi- 
cient functions are only polynomial in the quark masses. Only this fact allows for 
a consistent summation of logarithmic corrections through the choice /i^ = —q^. 
This task is accomplished in the next section. 
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However, the cancellation of mass log's only takes place up to the dimension 
of operators, to which all contributions have been included consistently. For 
example, the coefficient functions Cpp{q'^) still contain mass log's of the form 
Mm? In m, which were only cancelled, if mixing with operators of dimension 8 
had been included [|12|. 

The appropriate treatment of small masses in the OPE therefore is to expand 
in the mass up to operators of the dimension in question. Analogously, heavy 
masses have to be expanded in 1/M. The case of heavy quark masses will be 
discussed in a subsequent publication p7 . 



8 Renormalization Group Improved Coefficients 



The coefficient functions for the unit operator of eqs. ( |7.1| ) and ( [7.2| ) satisfy an 
inhomogeneous renormalization group equation. 



d 



where we have set Q"^ = — q^, and 



S,P 


N r 


— 


4vr2 [ 


V,A 


^ \ 


— 





- 3(M 



4 I 4\ 

m — 



(8.1) 

(8.2) 
(8.3) 



The presence of the inhomogeneity originates from the divergence of the current 
product of eq. ( p.l| ) 

The solution to this equation is given by 



ci{Q'/^?) = ciii) + ^ 

Po 



This expression is the improved form for the coefficient function C(, where the 
leading ln(Q^//i^) have been summed up. Po is the leading coefficient in the 
expansion of the /5-function, 



a.. 



/3(a,) = (3o— + 



Po = -^(lliV-2/), (8.5) 
5 



dn n 
where / denotes the number of quark flavours. Here and in the following we have 



adopted the notation of Pascual and Tarrach [^] for the renormalization group 
functions. 
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To improve on the contribution from higher dimensional operators, we note, 
that the sum in eq. ( p.ll ) has to be independent of /i. Therefore, the choice /i^ = 

sums up the leading log's, leaving us with condensates {Oi{Q^)), evaluated 
at the scale Q^. These can now be expressed in terms of condensates {Oi{fil)), 
evaluated at some fixed scale /xq, at which the numerical values of the condensates 
are known, e.g. /xq = 1 GeV. 

For simplicity, in the following we consider only one quark flavour q, however 
keeping the / dependence in the renormalization group functions. The general- 
ization to an arbitrary number of flavours should be obvious. 

In order to write the relation between (Oj(Q^)) and (Oj(/iQ)), it is convenient 
to assemble the operators Oj into a vector O, 



^ IT 



m 



(8.6) 



in which we have included the mass as an operator. This allows us to write the 
scale invariant combinations of the operators Oi in the form 



■>i, vo 



i?(/i) (0(/.)) . 



00 just is the invariant quark mass, the invariants 0i and 02 have been calculated 
in ref. , where also next-to- leading order corrections can be found, and 03 was 
obtained in ref. without taking into account mixing with the mass operator. 
The matrix R{fi) is given by 









V 












a. 



where 7^, 7°, and 7° are the leading order anomalous dimensions of the mass 
operator, the vacuum energy, and the mixed condensate, respectively. 



In 



3(Ar2 



70° 



N 
2^' 



and 



la 



{N^ - 5) 



.9) 



AN ' 27r' AN ' 

and (i° are defined to be the ratios = 7^//5o and d^ = la/Po: and finally, 
w, X, y, and z are given by 
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w 



70° ^ 3iV^ 



X = Z!^ = (8 10) 

(/?o + 7[L-7°) (8iV2 - 2iV/ - 3) ' ^ • ^ 

47^ (/5o + 7^ - 7° - 1) -6(iV2 - l)[8iV2 + 2iV(3 - /) - 3] 



(7° - I'JiPo + 7°. - 7°) (N^ + l)(8iV2 - 2iV/ - 3) 

-18A^2(jY2_i)[8jY2_^2A^(3-/)-3; 



(;3o + - 7^) 7r(iV2 + 1) _ 2iV/ - 3) (lOA^^ + 2Nf - 15) ' 

The element (1,4) of the matrix R which describes the mixing between [gsqcrFq) 
and is new. We should remark that strictly speaking it is not fully consistent to 
take the complete leading order matrix R at the order we have given the coefficient 
functions, because some entries of R appear first at 0{as)- Nevertheless we 
found it convenient to have at hand the full leading order renormalization group 
invariant combinations of the operators Oj. Use of the leading as well as next- 



to- leading order corrections will be made in |]T6|. 

The relation between condensates evaluated at two different scales is then 
given by 

(d(g^)) = uiQ',(,l){di(,l)), with f/(g^/i^) = r-\q') R{fj.l) . (8.11) 

The dependence on m{Q'^) appearing in the evolution matrix U^Q"^, /ig) can of 
course be expressed in terms of m{fiQ), using the relation 

such that the dependence of the evolution matrix enters only through ^^(Q^). 

Putting everything together, our final result for the renormalization group 
improved coefficient function of eq. ( |2.1| ) takes the form 



Po 



+ (7^V = Q')f>(Q',/"o)(0(/ig)), (8.13) 
where the appropriate powers of l/q^ have been absorbed into . 
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9 Summary 



We have calculated the coefficient functions of the quark and the mixed quark 
gluon condensate for scalar, pseudoscalar, vector, and axialvector current corre- 
lators to all orders in the quark masses, in the framework of the operator product 
expansion. For completeness the coefficient functions for the unit operator as 
well as the gluon condensate are reviewed. 

The proper factorization of short and long distance contributions has been 
performed, which requires renormalization of the condensates to the order consid- 
ered. It is found that it is only possible to absorb all long distance contributions 
into the condensates, if the OPE is expressed in terms of non-normal-ordered 
operators. This fact necessitates the inclusion of mixing with the unit operator 
under renormalization. 

The resulting coefficient functions were improved with the help of the renor- 
malization group equation. This is trivial in the case of the unit operator, but is 
somewhat more complicated for higher dimensional operators, since mixing under 
renormalization has to be taken into account. The scale invariant combinations 
up to operators of dimension 5, which are needed for the renormalization group 
improvement of the condensates, are given. The inclusion of mixing with the unit 
operator in the scale invariant of dimension 5 is new. 

An application of our results to an improved determination of the current 
strange quark mass, as well as a discussion of higher order corrections will be 
presented in a forthcoming publication |T6 
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Appendices 



A Expansions in one small mass m 



Here we list the expansions of the coefficient functions Cj in m, up to operators 
of dimension 6. Since the result in the case of Ci is very lengthy, it will only be 
presented up to dimension 4. With W — — q^, we find for the 
Scalar and Pseudoscalar currents: 

^4 ]\J2 ,,2 
-,2~ 



S,P 



N 



87r2 ^ 
^2mM 



2g^-3M^ + ^ln^-(2M^-g^)ln| 



In h In — 

o2 w W 



- 2m^ 


[l + ln^l 




. W. 



±2- 



w ' 

q2 



2 

, m 



+ 



m 



4 r 



o In 

.2 ^ 



(2M2 - g2) In 



m 



2^ 



(A.i: 



S,P 
QQ 

S,P 



Mq^ mq^i2M^-q^) , w?M^q^ m^M^q^ 4,,,^, 
± H ^7777^ ± 777; \ 7771 V 0[m )(A.2) 



W ' 2W^ 
— If m ± — + C>(m^) , 



1^3 

\q^\ > 



(A.3) 



12W\ 



f Mg2 



2W^ ' MW^ 



q^ + + In 



1^ 



S,P 
qFq 

S,P 
QFQ 



2 4 

m q 

Mq^ 
2 



g2 ^ 7^2 ^ g^2 1^ . 



mM 



W 



2 Jo 



2W^ 

-,2\ 



+ ©(m") 



Vector and Axialvector currents: 



N \ 10 



247r2 \ 3 



(A.4) 

(A.5) 
(A.6) 



— + 4M^ - 4— + 2(2M^ - 3g^) — In — + 2g^ In ^ 



1 + 2^-2^1^ 



-3 



r(2M2-g2)2 



c 



q- q- W \ "W^ 
-2(2M^-3a^ln — + 2gMn- 
mq^ 2m3g4(4M2 - g2^ 



+ 



31^4 



W 

(A.7) 
(A.8) 
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- M-'-$^'-f^M), l/l>M' (A,9) 



+ 0{m^) (A. 10) 



B The derivation of eq. (6.4) 

We consider the non-local vacuum expectation value {'■ gs(la{0)F^u{0) Ipi^x) :), from 
which we want to project out the local mixed condensate (:5's^a(0)F^i/(0) g/3(0) :). 
Taylor expansion of the quark field yields 

oo 
fc=0 

(B.l) 

where we have used the properties of the fixed-point gauge to express the ordinary 
derivatives in terms of covariant derivatives (see |^ and refs. therein). However, 



since we are not interested in the contributions from higher dimensional operators, 
we are free to replace the covariant derivatives by the ordinary ones, and work 
with the free equation of motion {i ^ — m)q{x) = 0. 

An Ansatz for (: (7sgQ,(0)F^i.(0) d^^ . . . d^^qpiGi) :) has to satisfy two conditions 
which come from the Lorentz-structure of eq. ( p.l| ): it has to be antisymmetric 
in /i and z/, but totally symmetric in /ii . . . /i^, due to contraction with the totally 
symmetric tensor x^^ . . . . 



Let us split the sum in (p.l| ) into "even" and "odd" terms, in order to simplify 
our notation. 



oo 



n=0 K'^'"^)- 

oo 

„=1 -LJ- 

where we have only explicitly kept the x-dependence. 

Consider first the even terms, k = 2n. There are two possibilities to construct 
Lorentz-invariant tensors which fulfill the above requirements. Therefore, our 
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Ansatz has the following general form 

{■■gsqaF^u •• • ^M2n 9/3(0) :) = A2n{(rf,u)i3aS,.i...,.2n (B.3) 

2n 

+ -B2n ^ I ('7^i/iJ^cj'S'^l.../ii.../i2nl' "I" i'^t^iu) I3a^fllli...fii...ll2n } • 

i=l 

Here 5'^i...^2n denotes the totally symmetric tensor which consists of a sum of A^2n 
products of n metric tensors, e.g., 

A hat on an index means, that this index is to be omitted. 
By induction it can be seen, that S^^,,,^^^ contains 

7V.„ = <M (B.5) 



elements ||3^. This number is important for the substitution of ( [B.3| ) into 
because it appears upon contraction of the two totally symmetric tensors >S'^i...^2n 
and x^^ . . . x'^^" , 

= Ar2„x2". (B.6) 

Contracting eq. (|K2|) with {(Tf"')^pS^'^-''^" and {a^'^'^)a|3S''^''■■■^''" , and using 
the equation of motion, leads to the system of equations 

DA2n + ^nB2n = r(D/2) (-m^)" i^.g^qaFq-) , 

4(D - l)2"r(n + D/2) ^ '^/' 

A2n + {D + 2n-2)B2n = 0, (B.7) 
whose solution is 

A - P + 2n-2)r(D/2)(-m^r F..-^ 
8(D - 1)(D - 2)2"r(n + 1 + D/2) ' 

8(L> - 1)(D - 2)2"r(n + 1 + D/2) ^ - ^ ' 

Bo = 0. 

Let us now come to the odd terms of eq. ( [B.2| ), k = 2n — l. Here the symmetry 
considerations allow for the following Ansatz, 

(B.9) 

+ -^271-1 I {,1 nS ^^...^2n-iy) f3a ~ iliy^ lit^h..t^2n-i) I3a } ' 
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The tensor S'^i...^2n-i is the sum of N2n-i products of one 7-matrix with n — 1 
metric tensors, where 

{2n-l)\ 
fn- 11 

As an example, for n = 2 



^-1 = ^^^^Tj:-^- (B-io) 



Similar contractions as in the even case lead to 

i TiD l2){—m'^Y 

A2n-l-lB2n-l = 0, (B.12) 



so that 



ir(D/2)(-m2)" 
~ Am{D - 1){D - 2)2«r(n + D/2) ^'^^^^^^-^ ' 



Putting now things together, i.e., substituting (IRSD , ( TO , (pIOD , and (pl3| ) 



into ( [B.2| ), and carrying out the contractions of x^^. . . x^" with the corresponding 
symmetric tensors, we arrive after some algebra at the simple result of eq. 
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